FLOW VISUALIZATION STUDIES 
OF THE XFV-12A 



Charles Lewis Peterson 



JUL 



Mi S..HOU1 

♦ .lA ^ 






























Monterey, California 







Approved for public release; distribution unlimited. 



T 167564 



UNCLASSIFIED 



SECURITY CLASSIFICATION of THIS page Date Entered) 



\ REPORT DOCUMENTATION PAGE 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


*. REPORT NUMBER 


2. GOVT ACCESSION NO. 


3. RECIPIENT'S CATALOG NUM8F.R 


4. TITLE (and Subtitle) 

Flow Visualization Studies of the 
XFV-12A 


S. TYPE OF REPORT ft PERIOD COVERED 

Master's Thesis; 

March 1975 


6. PERFORMING ORG. REPORT NUMBER 


7. AUTMORf*; 

, Charles Lewis Peterson 


«. contract or grant numbers; 


[S. PERFORMING ORGANIZATION NAME AND ADDRESS 

r Naval Postgraduate School 
Monterey, California 93940 


10. PROGRAM ELEMENT. PROJECT, TASK 
AREA A WORK UNIT NUMBERS 


11. CONTROLLING OFFICE NAME AND ADDRESS 

Naval Postgraduate School 
Monterey, California 93940 


12. REPORT DATE 

March 1975 


13. NUMBER OF PAGES 

92 


14. MONITORING AGENCY NAME a ADDRESS (U dl iterant from Controlling Office) 

Naval Postgraduate School 
Monterey, California 93940 


15. SECURITY CLASS, (of thie report) 

Unclassified 


15*. DECLASSIFICATION/DOWNGRADING * 

SCHEDULE ) 


16. DISTRIBUTION STATEMENT (of thie Report) 

^Approved for public release; distribution unlimited. 

j 


17. DISTRIBUTION STATEMENT (of the abstract entarad in Stock 20, if different from Raport) 


16. SUPPLEMENTARY NOTES 


13. KEY WORDS (Con'/nu* on revaree aide if nacaaaary and identify by block number) 

Recirculation 

Reingestion 

XFV-12A 

Vortex Shedding 


20. ABSTRACT (Continue on reverae aide if neceaaary and identity by block number) 

The XFV-12A incorporates thrust augmentation to generate 
propulsive lift. Under the influence of forward flight the 
lift jets are subject to bending, which produces shed vor- 
tices that induce velocities at the aircraft. Additionally, 
operations in ground effect during hover and for small for- 
ward velocities can result in reingestion of engine exhaust 
gases. The objectives of this thesis were twofold. The 



DD FORM 

w \ JAN 73 

(Page 1) 



1473 



EDITION OF I NOV 65 IS OBSOLETE 
S/N 0102*014* 660 1 | ^ 



UNCLASSIFIED 



SECURITY CLASSIFICATION OF THIS PAGE Dmt a Entared) 



UNCLASSIFIED 

CiiCUKiTY Classification of this D*r* 



Block #20 continued 



first was to investigate analytically aspects of the vortex 
shedding phenomenon through the application of potential 
flow around the lifting surfaces of the aircraft. Secondly, 
experimental analysis was performed to determine the unique 
reingestion properties of the XFV-12A operating statically 
under various wing and canard configurations. 

The pressure coefficient for jet curvature obtained 
under two-dimensional potential flow was 8.66. Reingestion 
of exhaust gases was primarily through . the main inlet and 
not the auxiliary inlet; the source of most of the exhaust 
gas was the canard. 



J 



DD Form 1473 (BACK) 
, 1 Jan 73 

S/N 0102-014-6601 



2 



SECURITY CLASSIFICATION OF THIS F* AGE(TFh».n Dt(» £nt»r»d) 



UNCLASSIFIED 



Flow Visualization Studies 
of the XFV-1ZA 



by 



Charles Lewis ^Peterson 

Second Lieutenant, United States Marine Corps 
B.S.A.E., United States Naval Academy, 1973 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN AERONAUTICAL ENGINEERING 



from the 

NAVAL POSTGRADUATE SCHOOL 
March 1975 



* 

V« j 

C * 



ABSTRACT 



The XFV-12A incorporates thrust augmentation to generate 
propulsive lift. Under the influence of forward flight the 
lift jets are subject to bending, which produces shed vor- 
tices that induce velocities at the aircraft. Additionally, 
operations in ground effect during hover and for small for- 
ward velocities can result in reingestion of engine exhaust 
gases. The objectives of this thesis were twofold. The 
first was to investigate analytically aspects of the vortex 
shedding phenomenon through the application of potential flow 
around the lifting surfaces of the aircraft. Secondly, ex- 
perimental analysis was performed to determine the unique 
reingestion properties of the XFV-12A operating statically 
under various wing and canard configurations. 

The pressure coefficient for jet curvature obtained under 
two-dimensional potential flow was 8.66. Reingestion of ex- 
haust gases was primarily through the main inlet and not the 
auxiliary inlet; the source of most of the exhaust gas was 
the canard. 
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I. INTRODUCTION 



The XFV-12A is a V/STOL aircraft currently under develop- 
ment by the Columbus Aircraft Division of North American 
Rockwell. Its design incorporates high performance super- 
sonic operation in a fighter or attack mission as well as 
vertical and short take-off operations from proposed ships 
of the 1980's. A unique property of the XFV-12A is incor- 
poration of thrust augmentation when operating under propul- 
sive lift. Combining the effects of hypermixing nozzles and 
Coanda blowing, this aircraft is designed to take off ver- 
tically with a thrust- to-weight ratio of less than one. 

When it becomes operational, the XFV-12A will offer the first 
employment of thrust augmentation in military or civilian 
aircraft . 

Problem areas under investigation include the shedding 
of vortices from the exhaust jets and an experimental analysis 
of the XFV-12A’s recirculation properties. The vortex shed- 
ding phenomenon occurs when a downflow of the lift jets is 
modified by forward velocity, ground effect, or both. The 
resultant vortices are shed as a result of a bending of the 
lift jets as shown hypothetically in Figure 1. This bending 
is derived from a high pressure area underneath and a low 
pressure area above the jet. These vortices are generated 
in much the same manner as wing tip vortices produced in 
finite airfoil theory. The nature and influence of these 
vortices are the unknowns under investigation. 
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A potential flow solution to a propulsive lift system 
was performed by Albers and Potter in 1971 [Ref. 1]. This 
analysis assumed an exhaust profile in accordance with 
Spence's theory [Ref. 2] for a jet augmented flap. The 
McDonnell Douglas Corporation [Ref. 3] has conducted work 
on a jet-wing lifting surface theory using elementary vor- 
tex distributions. A solution for the exhaust profile it- 
self, however, is not known. 

The vortex shedding is a consequence of the pressure 
difference on the two sides of a bent jet. To gain insight 
into the vortex problem, a potential flow model was formu- 
lated to determine jet bending. The potential flow was 
solved by two methods. A computer program was used to solve 
Laplace's equation for two-dimensional potential flow. For 
this case, only the canard's influence was considered. A 
continuous solution was obtained using an electrostatic 
analogy. Conducting paper provided the flow domain. Laplace's 
equation was again applied, but under the influence of both 
the wing and canard in a two-dimensional analysis. 

Recirculation experiments provide valuable data which 
adequately depict the reingestion properties of the XFV-12A. 
Since this characteristic is not readily calculated, results 
are based purely on experimental work performed and the cur- 
rent literature available on the subject. In subsequent 
discussions, the terms reingestion and recirculation will 
be used to describe the same phenomenon. 

Considerable work has been performed in the area of 
recirculation properties of lift engines. Investigations 



of recirculation for lift engines of various configurations 
applied to various airframe shapes have been performed 
[Refs. 4,5]. An analysis to determine the applicability of 
small scale testing has been done [Ref. 6] . Reduction of 
reingestion was examined by applying several concepts to 
prevent recirculated gases from entering lift engine inlets 
[Ref. 7] . Conclusions have not always been consistent when 
analyzing the most significant contribution to recirculation 
or the optimal method of prevention. However, it was commonly 
held that reingestion is a unique property for each airframe 
and accompanying engine configuration. Thus the particular 
properties of the XFV-12A operating in the vertical or short 
take-off mode offer a unique problem. 

The reingestion problem was investigated experimentally. 

The literature was helpful in formulating and designing ex- 
perimental apparatus and procedures. For carrying out the 
experiments, a one-tenth scale semi-span model of the XFV-12A 
was used. In order to simulate a zero forward velocity, 
testing was performed statically. Recirculation properties 
of interest included results for both the STOL and VTOL 
configurations at varying heights above the ground. In all 
cases, the reingestion at both the main and auxiliary inlets 
was determined. 

The computer program and its background are discussed 
in Section II of this thesis. Section III contains the work 
done with the electrostatic analogy. The next section describes 
the arrangement of the experimental apparatus including supple- 
mental hardware necessary to supply appropriate simulation of 
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the lift engines. Section V covers experimental work per- 
formed preliminary to the actual reingestion tests. The 
reingestion experiments are described in Section VI. Ex- 
perimental results, conclusions, and recommendations follow. 
Graphical test results and test data are recorded subsequent 
to the text. 
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II. COMPUTER SOLUTION TO THE POTENTIAL FLOW EQUATIONS 
FOR TWO-DIMENSIONAL THRUST AUGMENTOR CANARD 



A. ASSUMPTIONS 

In formulating the problem, the mass flow rate of the 
air leaving the diffuser exit plane was assumed to be equal 
to the secondary air entering the diffuser. This assumption 
was essential for satisfying continuity in the control volume 
which encircled the canard. (Using an experimentally ob- 
tained augmentation ratio of 6.84 [Ref. 8], the actual ratio 
of exhaust at the canard exit to the secondary flow is about 
eight to seven.) 

The working fluid was assumed to be incompressible, ir- 
rotational, and two-dimensional. With the irrotationality 
assumption, Euler's equations of motion for fluids become 
homogeneous. When a single dependent variable representing 
the potential function is introduced, a single partial dif- 
ferential equation results. 



l!l + iii = n 

9x2 9 y 2 u * 



( 1 ) 



The conjugate function to the potential function, or the 
stream function, is derived from the continuity equation for 
incompressible flow. 



The equation for the 
puter program. 



iii + ifi = n 

9x2 + 3 V 2 u - 



stream function 



was 



( 2 ) 



solved in the com- 
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B. NATURE OF PROBLEM 



The equation for the stream function is elliptical in 
nature yielding a boundary value problem. The boundary 
conditions were of the Dirichlet type. Values of the stream 
function were specified on all edges of the domain. Values 
of the stream function were also assigned to the canard 
flaps [Ref. 9] . The boundary conditions were held constant 
on the upper and lower bounds. The rate of change of stream 
function values on the upstream and downstream boundaries 
was held constant. These boundary conditions duplicated 
freestream conditions. 

For simplicity, only the canard was considered in this 
analysis. It was positioned nine chord lengths from the 
front boundary of the domain to its leading edge. The rear 
of the domain was positioned about 40 chord lengths aft of 
the canard. At this position it was anticipated that the 
flow field disturbances generated by the canard would have 
returned to freestream conditions at the boundaries. The 
discretized domain is illustrated in Figure 2. 

C. DISCRETIZATION 

Since the canard is small compared to the size of the 
domain, a variable mesh was necessary to represent optimally 
the flow while economizing on computer storage requirements. 
It was anticipated that the direction of the streamlines 
would be changing most rapidly in the vicinity of the canard. 
The highest concentration of solution points was then near 
the canai'd. The coarsest distribution of solution points 
was located near the corner most distant from the canard. 
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The final distribution of solution points within the domain 
resulted in 201 distinct points. 

Using a Taylor series expansion, the differential equa- 
tion can be reduced to a difference equation as illustrated 
in Appendix A. Each solution point is expressed in terms 
of the four points adjacent to it. The resultant error on 
transition to difference equations is of the fourth order. 
Deriving a difference equation for each solution point re- 
sults in 201 simultaneous, linear, algebraic equations in 
201 unknowns. In matrix form, the problem is easily stated 
as 

A*x = B (3) 

where A is the matrix of coefficients for the unknown values 
of the stream function, x is the vector of unknown solution 
points, and B is a vector containing the boundary conditions 
assigned to both the perimeter of the domain and the canard 
flaps . 

Since most of the elements in A are zero, it became more 
efficient to utilize a subroutine requiring only the storage 
of the nonzero elements and as few zero elements as possible. 
Such a library subroutine in the IBM 360 system at the Naval 
Postgraduate School is "GELB." This subroutine employs 
Gaussian elimination techniques in solving a banded matrix. 
This reduced the storage requirements for this particular 
problem by a factor of about ten. Not only did "GELB" re- 
duce the time required to obtain a solution, but it provided 
for a finer mesh size within the domain for a given quantity 
of core space. 
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III. ELECTROSTATIC ANALOGY SOLUTION TO THE POTENTIAL FLOW 
EQUATION FOR A TWO-DIMENSIONAL THRUST AUGMENTOR WING 
AND CANARD SYSTEM 

A. ELECTROSTATIC ANALOGY 

The differential form of Gauss’ flux theorem as it 
governs electric fields states that the divergence of the 
field is proportional to the charge density [Ref. 10]. In 
vector notation: 

V-E = p/e Q . (4) 

Since the electric field, E, at any point is equal to the 
negative of the gradient of the electric potential, V, 

E = - VV. (5) 



Gauss' flux theorem can be rewritten as Poisson's equation. 

V 2 V = -P/e Q . (6) 

If Poisson's equation is assumed to apply in a region 
of zero charge density, the expanded two-dimensional electric 
potential equation can be written: 



8 2 V + 9^V = Q 



(7) 



3x 2 8y 2 

This equation is Laplace's equation as it applies to elec- 
tric fields [Ref. 10] . 

If the voltage is allowed to be proportional to the 
stream function, a plot of equipotent ial lines on some con- 
ducting medium will yield a continuous solution for the 
stream function [Ref. 11]. Additionally, the zero charge 
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density assumption is appropriate since it is analogous to 
the irrotationality requirement of potential flow. 

Similarly, the current flowing in the domain is 
analogous to the fluid potential function. Since circula- 
tion is related to the potential function by 

r = / d<f> (8) 

the value of the current can be used to measure circulation 
[Ref. 11] . 

B. EXPERIMENTAL APPARATUS 

The hardware necessary to provide such a solution in- 
cluded conducting paper of a known resistance, an electro- 
meter for use in measuring voltages and currents, a ten-volt 
d-c power source, and seven potentiometers (Figure 3). The 
potentiometers took a predesignated fraction of the ten-volt 
source and applied it to a particular component of the wing- 
canard system. A special pen which contained a conducting 
paint was used to draw the wing and canard onto the conduct- 
ing paper. 

The resistance of the paper was determined using a known 
voltage source and the electrometer which measured current. 
The electrometer was placed in series with the conducting 
paper and the voltage source. The resistance was measured 
to be approximately 10,000 ohms per centimeter. The poten- 
tiometers were selected so that their resistances would be 
much less than that of the paper, but large enough to keep 
a low amount of current flowing from the ten- volt source. 
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If the maximum allowable current in the battery was con- 
strained to be ten milliamps, then the potentiometers had 
to be on the order of 1,000 ohms. This satisfied both the 
current and resistance constraints. 

The voltage induced on the wing and canard flaps were 
assigned so as to provide for equal circulation around the 
forward and aft flaps of the wing and canard. In addition, 
voltages were constrained to provide a flow velocity in the 
canard and wing augmentors which was twice that of freestream 
conditions . 

The equal circulation requirement was implemented by 
providing equal but opposite current flowing in the two 
Coanda flaps. The flow velocity constraint was imposed by 
comparing the distance between the Coanda flaps on the paper 
with the distance between the streamlines upstream 'which 
stagnated on those Coanda flaps. 

The center ejector flaps were constrained to have zero 
current and a voltage value midway between those assigned 
to the Coanda flaps. These requirements resulted in a 
negligible lift contribution made by the center ejectors as 
well as an equal distribution of wing and canard augmentor 
inlet flow on each side of the center ejector. 

In conducting the experiment, however, the equai-but- 
opposite circulation requirement could be imposed on the 
canard only. The results shown in Figure 4 indicate the 
streamlines having the same voltage value as those induced 
on the wing and canard flaps. The ground is the only por- 
tion of the boundary depicted. 
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IV. EXPERIMENTAL DESIGN AND MODEL MODIFICATION 



A. PLUMBING AND FIXTURES 

Compressed air was supplied by two compressors rated at 
17.2 and 49.5 cubic feet per minute, respectively. One was 
equipped with a 100 gallon storage tank, while the other had 
an SO gallon tank. These compressors were installed in paral 
lei with three additional storage tanks with a capacity of 
400 gallons each. The total storage capacity of the system 
was 1,380 gallons or 1S1.3 cubic feet. 

The compressors were limited by design to a pressure of 
137 lb/in 2 . The actual operating pressure of the system was 
constrained to be no more than 125 lb/in 2 . For the planned 
operating velocity at the ejectors of 136 ft/sec, substantial 
ly more than 15 minutes of experimental run time at the de- 
sired operating conditions was available. 

A single regulator and stop valve connected the storage 
tanks to the plenum, which supplied the ejector hypermixing 
nozzles and Coanda blowing. This provided for single-valve 
control of the experiments. Once the desired pressures had 
been obtained by altering the setting on the plenum valves, 
the plenum valve settings could be held constant. Only the 
stop valve was needed to initiate and terminate a particular 
experimental run. 

In anticipation of testing in the Naval Postgraduate 
School smoke tunnel, a two-inch pipe was channeled to the 
vicinity of the smoke tunnel from the storage tanks. In 
parallel with the outlet in the static testing area, the 
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pipe routed to the smoke tunnel was sealed during static 
testing. Simple transition to smoke tunnel testing can be 
performed by sealing the outlet in the static testing area 
and moving the regulator and stop valves to the smoke tun- 
nel outlet. 

The plenum has eight outlets of which only six are needed 
to supply the ejectors on the wing and canard. Three-quarter- 
inch-diameter automotive heater hoses were used to connect 
the plenum to the Coanda slots for both wing and canard. 
One-inch hoses were needed for both center ejectors (Fig. 

5). 

B. INLET SIMULATION 

In order to depict adequately aircraft operation in the 
static, position, some means of duplicating the effect of 
the inlets needed to be accomplished. Upon commencing work 
on the XFV-12A, only the main inlet had been fabricated 
with no means of providing suction for inlet flow. 

The auxiliary inlet was fashioned from copper sheet metal 
and attached to the top of the aircraft model as shown in 
Figure 6. For both the main and auxiliary inlets, ducts 
were constructed for the purpose of establishing the correct 
inlet flow rate. The ducts are shown in Figure 7. A suf- 
ficient quantity of inlet flow is pumped to satisfy the 
continuity equation as applied to inlet and exit conditions. 
Continuity required that the sum of the mass flow rates 
entering the main and auxiliary inlets was equal to the mass 
flow rate of primary air from the wing and canard. (For 
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purposes of these experiments, the contribution of fuel to 
the mass flow rates of the exhaust jets was assumed to be 
negligible . ) 

A first attempt at entraining ambient air through the 
inlets was done using the ejector principle. Copper tubing 
with six ejectors was fabricated to blow air along the main 
inlet channel nesr the face of the inlet. A photograph of 
the ejector appears in Figure 9. A similar copper ejector 
was constructed for the same purpose in the auxiliary inlet. 
These ejectors were found to function but could not provide 
sufficient mass flow to match that leaving the Coanda slots 
and ejectors of the propulsive system. Thus, some other 
means of duplicating the suction effect at the inlets had 
to be devised. 

A suction blower was attached to the inlet channels and 
sealed so that all of the entrained air was expelled out the 
exit plane of the blower; see Figure 10, Using this method 
the dynamic pressure measured at the exit plane of the 
blower was found to be 0.6 inches of water. For purposes of 
analysis of the continuity equation, the flow was assumed to 
be incompressible between the exit plane of the blower and 
the inlets. With this simplification, velocities and cor- 
responding dynamic pressures at the inlet and the blower 
exit plane were dependent only on the relative areas. 




(9) 



and 





( 10 ) 
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The dynamic pressure at the exit plane of the motor dic- 
tates the conditions which must exist at the ejectors for 
continuity to be satisfied for the system. 

Upon comparison of geometries, the required velocity at 
the main inlet was found to be 27.5 ft/sec. The correct 
velocity at the ejectors and Coanda slots of the wing and 
canard was then 136 ft/sec. Applying the augmentation 
ratio determined previously, the mean velocity at the exit 
plane of the diffusers must be 26.8 ft/sec. 

C. MODEL MOUNTING 

Once the inlets had been constructed, the model was 
mounted on three-quarter-inch plywood so that a plane of 
symmetry was established for testing purposes. The plane 
of symmetry is commonly known as a reflection plane. The 
other side of the reflection plane was cut axvay at the in- 
terior of the model to allow for the introduction of air 
from the plenum. A photograph of the reflection plane 
mounting at the static testing site is shown in Figure 8. 

The ground was simulated by another piece of plywood 
mounted perpendicular to the model plane of symmetry and 
parallel to the ground. This piece was attached to the 
reflection plane with bolts and could be moved up and down 
with respect to the aircraft at one-inch intervals. The 
upper limit was analogous to placing the aircraft on the 
ground with the gear lowered. The lower limit simulated 
the aircraft at an altitude of 20 feet and was accomplished 
by removing the ground plane entirely. 
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As the "ground' 1 was moved down from the aircraft during 
testing, it became necessary to use a larger piece of ply- 
wood to compensate for the enlarged area which was under 
the influence of the aircraft's exhaust jets. This was done 
by placing a larger piece of plywood covering 32 square feet 
directly over the original "ground." 
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V. PRELIMINARY EXPERIMENTAL WORK AND ANALYSIS 



A. NATURE OF RECIRCULATION 

Reingestion of hot exhaust gases is a problem of particu- 
lar interest to designers of VTOL aircraft and STOL aircraft 
when operating at forward speeds of less than 25 knots [Ref. 
12]. This phenomenon involves the entrainment of some por- 
tion of one or more exhaust jets into the intakes of the 
lift engines. Adverse effects on performance include a de- 
crease in the available thrust, large increase in inlet 
temperature, and distortion in the temperature profile at 
the inlets [Ref. 7]. These last two effects can lead to 
stalling of the compressor and engine surge [Ref. 7] . Ver- 
tical takeoff and landing operations over areas with loose 
terrain also pose a problem of ingestion of particles kicked 
up by the exhaust jets [Ref. 13] . Potential loss of avail- 
able control power was illustrated dramatically by a XC-142. 

A landing accident in 1965 was attributed to a loss in con- 
trol power due to recirculation of hot exhaust gases [Ref. 

12 ] . 

Exhaust gas reingestion has been classified according to 
the method in which exhaust gas is transported to the vicinity 
of the inlet. The first type is a diffusion of the exhaust 
gases into the ambient air. The buoyant force of the low 
density exhaust carries it up from the ground to the vicinity 
of the inlet. This is commonly called the "buoyant" effect 
[Ref. 5]. 
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The second type of recirculation involves an inter- 
action of the exhaust jets with the ground and subsequently 
with one another so that exhaust gases are forced upward in 
a "fountain" effect between the two exhaust jets. This type 
usually requires less time for its effects to be felt. The 
fountain effect usually is considered to be the most serious 
of the two types in terms of potential hazard to aircraft 
operation [Ref. 5]. 

The reingestion problem is one which cannot be treated 
analytically at present. Rather, experimental results are 
the only reliable source of information concerning the re- 
circulation characteristics of a particular airframe shape 
and engine configuration. Even the interpretation of re- 
sults has resulted in conflicting conclusions. For example, 
independent tests arrived at opposite conclusions with re- 
gard to the most significant type of reingestion transport 
[Refs. 5,7]. 

B. ANALYSIS OF VALIDITY OF SMALL-SCALE TESTING 

In order for experimental results to be valid, reinges- 
tion characteristics observed for the scale model must be 
applicable to the full-scale aircraft. Also, since operating 
conditions of the aircraft may differ from those imposed on 
the model during testing, reingestion must be evaluated as 
a function of differing exit temperature and pressure 
conditions . 

In an investigation of scaling effects on VTOL recircu- 
lation, recirculation was determined to be scalable both 
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statically and dynamically [Ref. 6] . (The experimental 
variables used were temperature, temperature gradient, and 
inlet distortion.) Another analysis determined reingestion 
to be independent of the temperature and pressure profiles 
exhibited by the engine over the range tested [Ref. 4] . 

These conclusions indicated that the experimental test con- 
ditions provided significant results even though full scale 
operating conditions were not duplicated. 

With specific regard to the XFV-12A, an evaluation of 
the characteristics of thrust augmentation indicated that 
the augmentation ratio of a system is almost independent of 
nozzle pressure ratio [Ref. 14] . Consequently, conditions 
unique to propulsive lift employing thrust augmentation were 
duplicated at the test conditions. 

C. SMOKE AND LASER TESTING 

The first experiments performed to determine the existence 
of recirculating exhaust gases were done using smoke intro- 
duced into the flow and a portable laser for detection. A 
smoke jet was produced from an oil smoke generator similar 
to those used in the NPS smoke tunnel. Using an air flow 
from one of the extra outlets of the operating plenum, the 
smoke was forced through a narrow length of plastic tubing 
to the vicinity of the canard ejector. 

The smoke was laminar for about 2 or 3 inches upon 
leaving the tube, followed by transition to turbulence as 
it was entrained with secondary air into the canard. Sub- 
sequent to entrainment the smoke A\ras impossible to trace 
without the use of the laser. 
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The smoke experiments were successful in demonstrating 
a small degree of recirculation, although no quantitative 
results could be determined. They also demonstrated that 
the main and auxiliary inlet ejectors functioned as desired 
but were unable to develop sufficient amounts of inlet flow 
to match the wing and canard ejectors. The results of smoke 
testing were inconclusive; however, the results did indicate 
that further investigation was desirable. 

D. EXPLORATORY FREON TESTING 

Experiments conducted in similar fashion to the smoke 
tests were performed using an aerosol can containing freon 
instead of smoke. A freon leak detector normally used to 
detect refrigerant leaks was used to measure reingestion. 

These results confirmed that exhaust gases were being 
re ingested, though there was again inadequate suction at 
both intakes. This experiment demonstrated the need for a 
better method of inlet simulation as well as further analysis 
of reingestion characteristics. 
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VI. REINGESTION EXPERIMENTS 



A. DESCRIPTION OF LEAK DETECTOR 

As a detection device, a General Electric refrigerant 
leak detector was acquired. Normally used to detect leaks 
as small as 1/2 oz/yr in coolant systems, it was aptly 
suited for use in detecting exhaust gases contaminated with 
freon. The detector is illustrated in Figure 11. 

The sensitivity of the detector was varied by a knob 
on the face of the instrument and by a switch which indi- 
cated either a high or low sensitivity setting. The probe 

i 

held a light which flashed when freon was absorbed. The 
flash rate gave a qualitative indication of the absorption 
rate, although the flash frequency was not proportional to 
the freon contamination of the mass flow being absorbed. 

A moderate rate of reingestion was defined as that 
flash frequency which was equal to that obtained when the 
probe was placed directly over the reference leak. Since 
the reference leak has a known mass flow rate, it can be 
used as a reference value. Three other ordinal quantitative 
descriptions, heavy, light, and trace, were used to repre- 
sent deviations from the reference "moderate" condition. 

Before each experimental run, the detector was adjusted 
by reducing or increasing the sensitivity until at free- 
stream conditions, the probe light just ceased to flash. 

Thus each set of data was acquired using the same reinges- 
tion descriptions. Contamination of the ambient air at low 



31 



levels of freon forced experiments to be performed at sub- 
stantial time intervals. This allowed the background 
freon contamination to dissipate and provided time to re- 
supply the compressed air storage tanks. 

B. EXPERIMENTAL DEPENDENT VARIABLES 

Of interest in performing the experiments was the flow 
path of recirculating gases from the exit plane of the dif- 
fuser of both the canard and wing to the main and auxiliary 
inlets. It was anticipated that these recirculation char- 
acteristics would vary with height above the ground and 
with the angle between the chord of the wing and canard 
flaps and the ground. The investigation covered recircula- 
tion characteristics as height above the ground varied from 
zero to about 24 inches (20 feet in full-scale dimensions). 
The angle of the canard and wing flaps was varied for each 
height at 10 degree increments from the vertical to an 
angle of 30 degrees from the vertical. Thus the pattern of 
reingestion was investigated for the aircraft operating in 
both the STOL and VTOL modes. 

In addition, a profile of the distribution of recircula- 
ted exhaust at the main inlet was determined. This was done 
by converting the six main inlet ejectros previously employed 
unsuccessfully to the role of gas sampling ports. Tests were 
performed at a height above ground of zero and with the wing 
and canard flaps in the VTOL configuration only. Profiles 
were acquired which reflected the distribution of exhaust 
gases from the front or rear Coanda slots or the ejector in 
the canard. 
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C. QUANTITATIVE REINGESTION ANALYSIS 



Some means were necessary of quantitatively specifying 
the meaning of the descriptions of degree of reingestion. 

Since the term "moderate" has been assigned to that amount 
of intake corresponding to the reference leak, it provided 
the only basis for a quantitative analysis. 

A result can be obtained analytically in a straightfor- 
ward manner if some assumptions are made concerning the con- 
ditions under which the experiments were conducted and the 
nature and properties of the freon. First, it was assumed 
that thermodynamic properties could be evaluated at 70°F. 

The freon was assumed to have the molecular weight of Freon- 12 
(CCI 2 F 2 ) ; see [Ref. 15] . The freon in the aerosol can was 
assumed to be a mixture of vapor and liquid at saturation 
conditions . 

Five mass flow rates must be calculated using known 
properties of freon and air and the previously outlined as- 
sumptions. The mass flow rate of freon from the reference 
leak is required; this is equal to -the mass flow rate for 
moderate recirculation. The mass flow rate of air in the 
diffuser was also needed, as well as the mass flow rate of 
the detection probe. Finally, the mass flow rates of the 
inlets and the aerosol can were calculated. The various 
mass flow rates and corresponding sources and sinks are 
diagramatically shown in Figure 12. 

The mass flow rate of the reference leak is known to be 
1/2 oz/yr from the instructions for use of the instrument. 

When the units are altered to yield lbm/sec, the result is 



33 



m Q = 0.99 x 10 9 lbm/sec. 

The mass flow rate of the aerosol can was 

the momentum conservation equation and solving 

Ij2 = 2 (AP) 

1 " <>f ' 



resolved us 
for veloci 

( 11 ) 



ing 

ty. 



The pressures and density were obtained from a table of 
properties of saturated freon [Ref. 15] . The calculated 
velocity is 



U x = 567 ft/sec. 

The diameter of the tube leading from the aerosol can was 
measured to be 0.025 inches. The mass flow rate was then 
readily calculated. 

m^ = 3.92 x 10 lbm/sec. 

The mass flow rate in the diffuser was computed using 
the results obtained for the augmentation ratio [Ref. 8] 
and the exit conditions of the wing and canard diffusers. 
Using a primary air velocity of 136 ft/sec, the mass flow 
rate of the canard diffuser is 

m 2 = 4.58 x 10 1 lbm/sec. 

The mass flow rate in the probe required a supplementary 
experiment. The probe was attached to a manometer at one 
end. The other end led to a reservoir of fluid whose upper 
surface was exposed to ambient conditions. For small changes 
in pressure, incompressibility can be assumed for the air. 

The flow rate of gas being drawn into the probe could be in- 
ferred by observing the rate of displacement of the manometer 
fluid. 
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Using the speed of the air in the manometer tube, the 
mass flow rate of the gas in the detector probe could be 
computed directly from the continuity equation for incom- 
pressible flow. Subsequent applications of these results 
assumed that air was the dominant gas being absorbed by the 
probe . 



U A = U, ,A. . . 
p p M M 



( 12 ) 



Solving for the only unknown, U , yields a velocity of gas 
in the probe of 5.37 in/sec. The mass flow rate of the 
probe was 

- e 



m 3 = 4.82 x 10 



lbm/sec. 



The mass flow rate of the main inlet was calculated 
using measurements taken at the exit plane of the electric 
blower. The dynamic pressure was recorded which yielded 
the velocity at the blower exit. Comparing geometries from 
the analysis described in Section IVB, the properties at 
both inlets could be calculated. The resultant mass flow 
rate for the main inlet was 



m^ = 6.29 x 10~ 2 lbm/sec. 

The mass flow rate of the auxiliary inlet was assumed 
to be the same as that of the main inlet. The area of the 
auxiliary inlet is greater than that of the main inlet. 
However, two ducts of equal area connected the inlets to 
the blower. These ducts were smaller in cross sectional 
than either inlet and consequently dictated the mass flow 
rate of each. 
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The first problem in quantitative analysis was the com- 
putation of the fraction of mass flow in the detection probe 
that was freon. For a moderate reingestion rate, this was 
represented by the ratio of the mass flow rate in the refer- 
ence leak to the mass flow rate of the probe. 

F 1 = m /m, = 2.06 x 10 4 . 

1 o 3 

The second problem was finding the concentration of 
freon in the diffuser of the canard. Since this concentra- 
tion was assumed to be constant from the point of introduc- 
tion to the ejector, it was found by computing the ratio of 
mass flow rate of freon introduced into the air hose to the 
mean flow rate of air at the exit plane of the diffuser. 

^2 = n! i/ m 2 = ^.58 x 10 3 . 

A reingestion percentage was then calculated yielding 
the concentration of freon being ingested as compared to 
that being expelled. For moderate reingestion conditions 

F = F 1 /F 2 = 2.40%. 

Consequently a moderate reading indicated that 2.40% of the 
diffuser exit mass flow was being reingested by the main 
inlet . 

Further, it was desirable to know the percentage of 
inlet flow that- was recirculated from the engine exhausts. 
This was computed simply by multiplying the reingestion 
percentage by the mass flow rate of the diffuser, then 
dividing by the mass flow rate of the main inlet. 

G = F x n^/m^ = 17.8%. 
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Therefore a moderate reingestion reading indicated 
that about 18% of the inlet flow was comprised of recircu- 
lated exhaust gases. Here exhaust gases were defined as 
the mixture of entrained secondary air and primary air ex- 
pelled from the canard or wing augmentor. 

D. EXPERIMENTAL PROCEDURE 

Freon was injected through drilled holes in the heater 
hoses leading to the Coanda slots and the ejectors. The 
aerosol can could be moved readily to introduce freon into 
any given flap. The six points of introduction were as- 
signed station numbers sequentially from the forwardmost 
canard flap aft. The flaps and their accompanying station 
numbers are shown in Figure 13. Holes were drilled in the 
ducts leading from the inlets to which the detector probe 
was attached. The wing and canard flaps were set by mea- 
suring the angle between the chord of each flap and the 
vertical. The stiffness of the hoses was sufficient to 
maintain the position of the wing and canard flaps through- 
out a set of experiments. 

The experiments were conducted by varying the height 
above the ground with a constant wing and canard angle. 
Measurements were taken at both the main and auxiliary in- 
lets for each of the six points of introduction. Thus 
twelve readings were taken at each height above the ground 
for each wing and canard angle. The wing and canard were 
rotated together, so that there were no experimental runs 
for which the wing and canard angles were different. A 
total of 288 reingestion readings were taken. 
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The results were plotted according to the twelve com- 
binations of inlet and point of introduction. The twelve 
reingestion profiles are shown in Figures 14-25. The 
abscissa of each graph is the angle between the chord of the 
wing and canard flaps and the ground. The ordinate is the 
height above the ground, non-dimensionalized with respect 
to a canard reference length. The canard reference length 
was defined as the square root of the area of the entrance 
plane of the canard augmentor. Its value is 4.90 in. 

The experiments covering the inlet reingestion profile 
were performed with the canard since the wing did not con- 
tribute substantially to main inlet recirculation. Five of 
the six different gas sampling ports in the inlet face were 
sealed with adhesive tape at any one time, so that the rein- 
gestion reading measured only streamtubes passing near the 
open port. The response of all six sampling ports on the 
inlet face was recorded for each of the three points of in- 
troduction. A total of 18 readings were taken. Patterns of 
the main inlet reingestion are shown in Figures 26-28. 
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VII. EXPERIMENTAL RESULTS AND INTERPRETATION 



A. POTENTIAL FLOW SOLUTION 

The digital potential flow solution yielded a result 
which indicated that "supercirculation" was present at the 
rearmost canard flap. Here supercirculation was defined to 
occur when a streamline circumscribed any particular lifting 
surface. Supercirculation was present throughout the range 
of variables of velocity in the diffuser and mass flow dis- 
tribution entering the canard augmentor. There was also a 
rapid return of the exhaust profile to freestream conditions 
only a few chord lengths aft of the canard. 

The electrostatic analogy solution eliminated inaccuracies 
in the coarse mesh size. The flow around the canard was ab- 
sent of supercirculation. However, a stagnation point was 
produced on the upper surface of the rear canard flap. The 
influence of the wing and its exhaust was apparent in its 
influence on the canard flow. The tendency of the canard 
exhaust jet to return to freestream conditions was negated 
by the presence of the wing and its flow field. The highest 
velocities in the flow field occurred between the wing and 
the simulated ground. The bending of the exhaust jets and 
resultant vortex shedding was evident primarily around the 
wing. Supercirculation was evident around the rearmost flap 
of the wing for this case. 

Using the results of the electrostatic analogy as an 
indication of the velocity and shape of the exhaust jets, 
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an analysis was done to determine the effect of the exhaust 
jets and vortex shedding on the lift of the aircraft. Using 
the results shown in Figure 4, the radii of curvature of the 
exhaust streamlines and the distances between the adjacent 
streamlines in the normal direction were measured graphically. 
This was done in the vicinity of the trailing edge of the rear 
wing flap. 

The rate of change of pressure in the normal direction 
is related to the local dynamic pressure and streamline radius 
of curvature by 



U 2 
'a Rij 



9JP 

9n 



(13) 



Considering streamlines one and two of Figure 4, the above 
equation can be integrated to yield the following result. 



P - P 
r 2 *1 



- P 



U' 



a R 



12 



(n. 



n-i ) • 



(14) 



Nondimensionalizing the pressure difference with respect to 
the freestream dynamic pressure, the following relation can 
be derived. 



(P 2 - P x ) _ 2(AnJ 2 

q co R 12 ^ n 2 " n l^ 



(15) 



The results of these measurements and a resultant pressure 
coefficient across the exhaust jet are shown in Appendix C. 

Further, since the lift coefficient can be represented 
by the difference in the average pressure at the top of a 
surface and the bottom of the surface divided by the dynamic 
pressure, the resultant value for AP/q^ across the jet is 
analogous to a lift coefficient. The numerical value for 
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the pressure coefficient across the jet is greater than that 
which would have been obtained from a three-dimensional ap- 
Pl ication of potential flow. 

B. RECIRCULATION 

Recirculation was most severe in the VTOL configuration 
when the aircraft was situated on the ground. The main in- 
let was the largest recipient of exhaust gases. In general, 
reingestion decreased as height above the ground increased. 
However, for some cases, re ingestion increased as height 
aobve the ground increased and then decreased again as the 
height was further increased; for example, see Figure 15, 

0 = 60°, and Figure 16, 0 = 90°. The highest readings re- 
sulted when the point of introduction was at the forward 
Coanda slot of the canard. 

The results of introducing freon into the canard or 
wing were not independent of the point of introduction. 
Usually there was considerable difference in the response 
from station to station within the wing or canard. For 
example, see Figures 14, 15, and 16 for the response of 
the main inlet to introduction points within the canard. 

The angle that the wing and canard made with the ground 
had considerable influence on recirculation. With the wing 
and canard at an angle of 60 degrees with the ground, the 
wing was never a source of reingestion. The canard pro- 
vided some reingested gases at all altitudes and configura- 
tions. Deviations in wing and canard angle with the 
vertical resulted in a substantial reduction of reingestion 
for all altitudes and configurations. 
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The distribution tests at the main inlet indicated that 
reingestion was concentrated most heavily on the side most 
distant from the plane of symmetry of the aircraft when the 
point of introduction was station 1; see Figure 13 for 
station identification and Figure 26 for the measured re- 
sult. For station 2, however, the reingestion properties 
were most pronounced at the upper portion of the main in- 
let; see Figure 27. When freon was introduced into station 
3, there was an even distribution over the entire main inlet 
face; see Figure 28. The level of reingestion was only 
trace for this case. 

The response time for all reingestion readings was 
small, indicating that the fountain effect was the only ef- 
fect contributing to recirculation. There was no opportunity 
for the "buoyant" effect to be felt, however, since there 
was no significant temperature difference between the pri- 
mary and ambient air. Responses occurring a considerable 
length of time after the introduction of freon were attri- 
buted to diffusion and mixing of freon into the ambient 
air. These responses were never greater in magnitude than 
trace . 
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VIII. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The potential flow solution was valuable in obtaining an 
approximation to the pressure coefficient across the exhaust 
jet. This coefficient provided insight into the contribu- 
tion of a two-dimensional propulsive lift jet to the flow 
field of the aircraft. 

Reingestion poses a problem primarily to the main inlet 
under VTOL conditions. The magnitude of the problem was 
greater than the 18 % obtained with moderate reingestion 
for all altitudes. The canard was the most significant 
source of reingested exhaust gases, and was the only source 
to supply readings above moderate. The wing contributed no 
gases for reingestion in configurations other than VTOL. 

B. RECOMMENDATIONS 

Further work would be desirable to determine the effect 
of a forward velocity as simulated in the wind tunnel on 
reingestion properties. Closer scrutiny might be applied 
to the flow path from the canard ejector and Coanda slots 
to the main inlet. 
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APPENDIX A: SAMPLE DERIVATION OF DIFFERENCE EQUATION 



Each linear, algebraic equation was derived as follows 

iKx,y+k) 

j} k 

*4»(x+j,y) 



iKx,y-k) 

beginning with Taylor series for each axis of the domain. 
t' (x,y+k) = \p (x ,y) + ^ (x,y) + — |yf- (x,y) + ... 

and 

^(x,y-k) = i|>(x,y) - (x,y) + |r fyr ( X >X) ' ••• 

Adding the above equations and solving for the second par- 
tial, the following equation results. 



= §r [ 4»(x,y+k) - 2iJ>(x,y) + iKx,y-k) ] + 0(4). 

Similar analysis for the second partial in the x-direction 
yields 

f^T = J T t ^C x+ 3 »y) ' 2tp( x ,y) + iKx-j,y) ] + 0(4). 



Applying Laplace's equation 

8 2 iJ> 



+ ifi = o 

8x2 8 y 2 u 



each solution point is represented as a linear function of 
its four adjacent points. 



<K x »y) = -nrk 



2' f j 2 + k ~ 2 ) ' [ Kx,y+k) + ^ (x,y-k) ] + 2 ( jV+kTJ 
[ ^(x+j ,y) + i|<(x-j ,y) ] . 
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APPENDIX B: MASS FLOW RATE CALCULATIONS OF QUANTITATIVE 

ANALYSIS 



The following data was obtained either experimentally 

or from tables of known properties of freon or air. 

P a = 7.66 x 10' 2 lbm/ft 3 

p f = 2.028 lbm/ft 3 

U-^ = 567 ft/sec 

^ = 136 ft/sec 

Up = 0.45 ft/sec 

U^ = 27.5 ft/sec 

A^ = 4.91 x 10 4 in 2 

A 2 = 0 . 80 7 in 2 

A = .0204 in 2 

P 

A^ = 4.3 in 2 



The mass flow rates used in the quantitative analysis 
were calculated using the above data and appropriate conver- 
sion factors to maintain consistency in units. 



m l = p f u l A l = 3,92 


X 


10 ' 3 


lbm/sec 


m 0 - p U 0 A„ = 4.58 
2 a 2 2 


X 


10' 1 


lbm/ sec 


m 7 =pUA =4.82 
3 ^a p p 


X 


10' 6 


lbm/sec 


m 4 = P a^4^4 = 29 


X 


10' 2 


lbm/sec 
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APPENDIX C: RESULTS OF POTENTIAL FLOW LIFT JET ANALYSIS 
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Figui'e 





. SHED VORTICITY ASSOCIATED WITH LIFTING WINGS 
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Resultant stream function values are shown in sample 
computer program output. 



Figure 2. DISCRETIZED DOMAIN. 
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Figure 3. ELECTROSTATIC ANALOGY EQUIPMENT. 
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Figure 5. MODEL SUPPORT, PLENUM, AND HEATER HOSES. 
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Figure 6. AUXILIARY INLET 
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Figure 7. REAR VIEW OF MAIN AND AUXILIARY INLET DUCTS 
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Figure 8. MODEL MOUNTED IN VTOL POSITION AT GROUND HEIGHT. 
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Figure 9. MAIN INLET EJECTOR. 
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Figure 10. BLOWER WITH ACCOMPANYING DUCTS AND DETECTOR PROBE IN PLACE. 
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Figure 11. LEAK DETECTOR, PROBE, AND FREON SOURCE. 




Figure 12. PATH OF RE INGESTED EXHAUST GASES. 
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(a) CANARD STATIONS 




(b) WING STATIONS 

Figure 13. IDENTIFICATION OF WING AND CANARD FLAPS. 
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Pigure 14. REINGESTION PATTERN; MAIN INLET, STATION 1. 
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Figure 15. REINGESTION PATTERN; MAIN INLET, STATION 2. 
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Figure 16. REINGESTION PATTERN; MAIN INLET, STATION 3. 
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Figure 17. REINGESTION PATTERN; MAIN INLET, STATION 4. 
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Figure 18. REINGESTION PATTERN; MAIN INLET, STATION 5. 
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Figure 19. REINGESTION PATTERN; MAIN INLET, STATION 6. 
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Figure 20. REINGESTION PATTERN; AUXILIARY INLET, STATION 1. 
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Figure 21. REINGESTION PATTERN; AUXILIARY INLET, STATION 2. 
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Figure 22. REINGESTION PATTERN; AUXILIARY INLET, STATION 3. 
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Figure 23. RE INGEST I ON PATTERN; AUXILIARY INLET, STATION 4. 
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Figure 24. REINGESTION PATTERN; AUXILIARY INLET, STATION S. 
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Figure 25. REINGESTION PATTERN; AUXILIARY INLET, STATION 6. 
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Figure 26. REINGESTION PATTERN; MAIN INLET FACE, STATION 1. 
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Figure 27. REINGESTION PATTERN; MAIN INLET FACE, STATION 2. 
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Figure 28. REINGESTION PATTERN; MAIN INLET FACE, STATION 3. 
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COMPUTER PROGRAM 



On the next page, the text of the computer program 
begins which solved Laplace's equation using finite dif- 
ference approximations. A total of nine different runs 
were made, each of which applied a different combination 
of velocity at the canard exit plane and ratio of mass flow 
on each side of the canard ejector. Each combination was 
obtained by altering the boundary values assigned to the 
canard flaps . 

The results of each of the runs were similar in shape. 
Consequently, only one of the outputs is included subsequent 
to the program. 
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ooo 



300 

400 

500 

900 



5 

10 



15 



20 



25 



30 



FORMAT (//) 

FORMAT (6F8.2»6F4.2»//> 

FORMAT (6F8.2,2F4.2,4X,3F4.2, / / ) 

FORMAT (/////) 

DIMENSI CN A (4947 ) , B < 20 1 , 1 ) , C ( 20 1 > 

DIMENSION D( 20 1 ) 

INPUT THE ELEMENTS OF THE COEFFICIENT MATRIX 

00 5 I = 1,4869 
A ( I ) = 0.0 
CONTINUE 
DO 10 I = 1,201 
C ( I ) = 0. 

CONTINUE 
M = 2 



= 13,69,14 



= 112,168,14 

N 

M 



N = 0 
DO 15 I 
M = M + N 
J = I + M 
K = J - 1 
L = J + 1 
A ( K ) = 1. 

A ( L > = 1. 

A ( J ) = -10. 

K = J + 12 
A ( K ) = 4. 

N = N -t- 1 
CONTINUE 
M = 2 
N = 6 
DO 20 I 
M = M + 

J = I + 

K = J - 1 
L = J + 1 
A { K ) = 1. 

A ( L ) = 1. 

A ( J ) = -4. 

K = J + 12 
A ( K ) = 1. 

N = N + 1 
CONTINUE 

DO 25 I = 223, 348,25 
‘ + 11 
+ 12 
+ 13 
+ 24 
‘ 4. 

= 1. 

: - 10 . 

> 1. 

■ 4. 

CONTINUE 

DO 30 I = 523,648,25 
' + 11 
+ 12 
+ 13 
+ 24 
= 4. 

= 1. 

= -10. 

: 1. 

; 4 . 

CONTINUE 

DO 35 I = 823,948,25 
11 
12 
13 
24 



J = I 
K = I 
L = I 
M = I 
A ( I ) i 
A ( J ) : 
A ( K ) 

A ( L ) 

A ( M ) 



J = I 
K = I 
L = I 
M = I 
A ( I ) 
A( J) 

A ( K ) 

A ( L ) 
A(M) 



J = I 
K = I 
L = I 
M = I 
A ( I ) = 4. 
A ( J ) = 1. 



+ 

+ 

+ 
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A ( K ) = -10. 

A < L ) = 1. 

A (M ) = 4. 

35 CONTINUE 

DO 40 I = 2248,2373, 
J = I + 11 
K = I + 12 
L = I + 13 
M = I + 24 
A ( I ) = 4. 

A ( J ) = 1. 

A ( K) = -10. 

A ( L ) = 1. 

AIM) = 4. 

40 CONTINUE 

DO 45 I = 3023,3123, 
J = I + 11 
K = I + 12 
L = I + 13 
M = I + 24 
A ( I ) = 1. 

A ( J ) = 4. 

AIK) = -10. 

AIL) = 4. 

AIM) = 1. 

45 CONTINUE 

DO 50 I = 3323,3423, 
J = I + 11 
K = I + 12 
L = I + 13 
M = I + 24 
All) = 1 . 

A ( J ) = 4. 

AIK) = -10. 

AIL) - 4. 

AIM) = 1. 

50 CONTINUE 

DO 55 I = 3623,3723, 
J = I * 11 
K = I + 12 
L = I + 13 
M = I + 24 
AID = 1. 

A I J ) = 4. 

AIK) = -10. 

AIL) = 4. 

AIM) = 1. 

55 CONTINUE 

DO 60 I = 3923,4023, 
J = I + 11 
K = I + 12 
L = I + 13 
M = I + 24 
All) = 1 . 

A ( J ) = 4. 

AIK) = -10. 

AIL) = 4. 

AIM) = 1. 

60 CONTINUE 

DO 65 I = 4223,4323, 
J = I + 11 
K = I + 12 
L = I + 13 
M = I + 24 




AIL) = 4. 

AIM) = 1. 

65 CONTINUE 

DO 70 I = 4523,4623, 
J = I + 11 



25 

25 

25 

25 

25 

25 

25 
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70 

75 

80 

85 

90 

95 



K = I + 12 
L = I + 13 
M = I + 24 
A ( I ) = 1. 

A I J ) = 4 . 

AIK) = -10. 

A ( L 1 = 4. 

A ( M ) = 1. 

CONTINUE 

DC 75 I = 398,498,25 
J = I + 11 
K = I + 12 
L = I + 13 
M = I + 24 
A < I I = 1. 

A ( J ) = 1. 

A ( K ) = -4. 

A ( L ) = 1. 

AIM) = 1. 

CONTINUE 

DO 80 I = 698,798,25 
J = I + 11 
K = I + 12 
L = I + 13 
M = I + 24 
A ( I ) = 1. 

A ( J ) = 1. 

A (K) = -4. 

AIL) = 1. 

AIM) = 1. 

CONTINUE 

DO 85 I = 2148,2223,25 
J = I + 11 
K = I + 12 
L = I + 13 
M ~ I + 24 
AID = 1. 

A ( J ) = 1. 

AIK) = -4. 

AIL) = 1. 

AIM) = 1. 

CONTINUE 

DO 90 I = 2848,2973,25 
J = I + 11 
K = 1 + 12 
L = I + 13 
M = I + 24 
All) = 1. 

A I J ) = 1. 

AIK) = -4. 

AIL) = 1. 

AIM) = 1. 

CONTINUE 

DC 95 I = 3148,3273,25 
J = I + 11 
K = I + 12 
L = I + 13 
M = I + 24 
AID = 1. 

A I J ) = 1. 

AIK) = -4. 

AIL) = 1. 

AIM) = 1. 

CONTINUE 

DO 100 I = 3448,3573,25 
J = I + 11 
K = I + 12 
L = I + 13 
M = I + 24 
AID » 1. 

A < J ) = 1. 

AIK) = -4. 
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A < L ) = 1. 

A ( M ) = 1. 

100 CONTINUE 

DO 105 I = 3746,3873,25 
J = I + 11 

K = I + 12 

L = I + 13 

M = I + 24 

A( I) = 1. 

A ( J ) = 1. 

A (K) = -4. 

A(U = 1. 

A ( M J = 1. 

105 CONTINUE 

DO 110 I = 4048,4173,25 
J = I + 11 

K = I + 12 

L = I + 13 

M = I + 24 

A < I ) = 1. 

A ( J ) = 1. 

A ( K J = -4. 

A(l) = 1. 



110 CONTI 

DO 115 I = 4348,4473,25 
J = I + 11 
K = I + 12 
L = I + 13 
M = I + 24 
A ( I ) = 1. 

A ( J ) = 1. 

A ( K ) = -4. 

Ad) = 1. 

ACM) = 1. 

115 CONTINUE 

DO 120 I = 2423,2523,25 
J = I + 11 
K * I + 12 
L = I + 13 
M = I + 24 
A ( I ) = 1. 

A ( J ) = 1. 

A ( K ) = -'4. 

ACL) = 1. 

ACM) = 1. 

120 CONTINUE 

DO 125 I = 1048,1098,25 
J = I + 11 
K = I + 12 
L = I + 13 
M = I + 23 
A ( I ) = 1. 

A ( J ) = 1. 

A ( K 5 = -4. 

AIL) = 1. 

AIM) = 1 . 

125 CONTINUE 

DO 130 I = 1123, 1248, 25 
J = I + 11 
K = I + 12 
L = 1 + 13 
M = I + 23 
AID = 6. 

A I J ) = 5. 

AIK) = -50. 

Ad) = 5. 

AIM) = 32. 

130 CONTINUE 

DO 135 I = 1398,1523,25 
J * I + 1 
K * I + 11 
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135 

140 

145 

150 

155 

160 



L = I + 12 
M = I + 13 
N = I + 23 
A ( J ) = 16. 

A ( K ) = 1. 

AIL) = -34. 

A ( M ) = 1. 

AIN) = 16. 

CONTINUE 

DO 140 I = 1598, 1648, 25 
J = I + 1 
K = I + 11 
L = I + 12 
M = I + 13 
N = I + 23 
A ( J ) = 4. 

A ( K ) = 1. 

A ( L ) = -10. 

AIM) = 1. 

A ( N ) = 4. 

CONTINUE 

A ( 166 1 ) = 0. 

DO 145 I = 1673,1798,25 
J = I + 1 
K = I + 11 
L = I + 12 
M = I + 13 
N = I + 23 
A ( J ) = 32. 

AIK.) = 5. 

A ( L ) = -50. 

A ( M ) = 5. 

A (N > = 8. 

CONTINUE 

DO 150 I = 1948,2073,25 
J = I + 1 
K = I + 11 
L = I + 12 
M = I + 13 
N = I + 24 
A ( J ) = 4. 

A { K ) = 1. 

A ( L ) = -10. 

ACM) = 1. 

A ( N ) = 4. 

CONTINUE 

DO 155 I = 2548,2673,25 
J = I + 11 
K = I + 12 
L = I + 13 
M = I + 24 
A ( I ) = 8. 

A ( J ) = 3. 

A ( K ) = -18. 

A ( L ) = 3. 

A ( M ) = 4. 

CONTINUE 

DO 160 1 = 2723, 2823,25 
J = I + 11 
K = I + 12 
L = I + 13 
M = I + 24 
A ( I ) = 2. 

A ( J ) = 3. 

A ( K ) = -9. 

A { L ) = 3. 

A(M) = 1. 

CONTINUE 
N = 5 
M = 0 

DO 170 I = 4648,4743, 19 
L = I + M 
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J = L + 11 
K = l + 12 
A(L) = 1. 

A ( J ) = 1. 

A ( K ) = -4. 

J = l + 13 
A ( J ) = 1. 

M = M + N 
N = N - 1 
170 CONTINUE 
N = 3 
M = 0 

DO 175 I = 4795,4837,14 
L = I + M 
J = L + 11 
K = L + 12 
A ( L ) = 1. 

A ( J ) = 4. 

A ( K ) = -10. 

J = L + 13 
A ( J } = 4. 

M = M + N 
N = N - 1 
175 CONTINUE 

A ( 1 ) = -10. 

A 12) = 1. 

A ( 13 ) = 4. 

A ( 99 ) = 1. 



A( 100) 
A( 101 ) 
A( 112) 
A ( 211 ) 

A ( 234 ) 

A ( 373 ) 

A (384) 

A (385 ) 

A ( 386 ) 

A ( 397 ) 

A ( 51 1 ) 
A( 534) 

A ( 673 ) 

A (684) 
A( 685) 

A (686) 

A ( 697 
A ( 8 11 ) 

A (834) 

A ( 973 ) 

A (984) 

A ( 985 ) 

A (986) 

A ( 997 ) 

A (998 ) 
A( 1009) 
A ( 1010 ) 
A ( 101 1 ) 
A( 1022) 
A ( 1023 ) 
A ( 1034 ) 
A( 1035) 
A( 1036) 
A ( 1 1 1 1 ) 
A ( 1 134 ) 
A ( 1298) 
A( 1309) 
A ( 1310 ) 
A( 1321 ) 
A( 1273) 
A( 1284) 
A( 1285 ) 
A( 1286) 
A( 1296) 



- 9 . 

2 . 

3 . 

0 . 

0 . 

3 . 

1 . . 
- 9 . 

2 . 

3 . 

0. 

0 . 

3 . 

1 . 

- 9 . 

2 . 

3 . 

0 . 

0 . 

3 . 

1 . 

- 9 . 

2 . 

3 . 

1 . 

= 1 . 

: - 4 . 

= 1 . 

= 1. 

• 4 . 

= 3 . 

= - 18 . 
= 3 . 

= 0. 

‘ 0. 

= 4 . 

= 3 . 

= - 18 . 
= 8 . 

= 8 . 

* 2 . 
- 30 . 

' 4 . 
16 . 
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A ( 1324 ) = 2. 

A ( 1335 ) = -12 
A (1336) = 2. 

A ( 1346 ) = 4. 

A (1349) = 4. 

A ( 1359 ) = 3. 

A ( 1360 ) = -18 
A ( 1361) = 3. 
A( 137 1 ) = 8. 

A ( 1374 ) = 4. 

A ( 1384 ) = 3. 

A ( 138 5 ) = -18 
A ( 1396 ) = 8. 

A ( 1409 ) = 0. 

A ( 1549) = 24. 
A ( 1559) = 2. 

A ( 1560 ) = -54 
A ( 1561 ) = 4. 

A ( 157 1 ) = 24. 
A ( 1574 ) = 3. 

A ( 1584) = 1. 

A ( 1585 ) = -9. 
A ( 1 596 ) = 3. 

A ( 1 609 ) = 0. 

A ( 166 1 ) = 0. 

A ( 1684 ) = 0. 

A ( 1 824 ) = 8. 

A ( 1834 ) = 4. 

A ( 1835 ) = -30 
A ( 1 836 ) = 2. 

A ( 1846 ) = 16. 
A ( 1 849 ) = 8. 

A ( 1859 ) = 3. 

A ( 1 860 ) = -18 
A ( 1871 ) = 4. 

A (1874) = 4. 

A ( 188 5 ) = -12 
A ( 1 886 ) = 2. 

A ( 1897 ) = 2. 

A ( 1899 ) = 8. 

A ( 1909 ) = 3. 

A ( 1910 ) = -18 
A ( 191 1 ) = 3. 

A (1922) = 4. 

A ( 1924 ) = 8. 

A ( 1934 ) = 3. 

A ( 193 5 ) = -18 
A (1947) = 4. 

A ( 1959 ) = 0. 

A ( 2 09 9 ) = 3. 

A ( 2 109 ) = 1. 

A ( 2110) = -9. 
A ( 211 1 ) = 2. 

A ( 2122 ) = 3. 

A ( 2 124 ) = 1. 

A ( 2 134 ) = 1. 

A ( 2 135 ) = -4. 
A ( 2 136 ) = 1. 

A ( 2147 ) = 1. 

A ( 2148 ) = 0. 

A ( 2236 ) = 0. 

A ( 2259 ) = 0. 

A ( 2398 ) = 3. 

A ( 2409 ) = 1. 

A ( 2410 ) = -9. 
A ( 241 1 ) = 2. 

A ( 2422 ) = 3. 

A ( 2536 ) = 0. 

A ( 2559 ) = 0. 

A ( 2698 ) = 4. 

A ( 2709 ) = 2. 

A ( 2710) = -12 



c 

c 

c 



A<2711) 
A ( 2722 ) 
A ( 2836 ) 
A ( 2859 ) 
A { 2998 ) 
A ( 3009 ) 
A ( 30 1 0 ) 
A < 301 1 ) 
A ( 3022 ) 
A (3136 ) 
A ( 3159 ) 
A ( 3298 ) 
A (3309) 
A (3310) 
A (331 1 ) 
A ( 3322 ) 
A ( 3436 ) 
A ( 3459 ) 
A (3598) 
A (3609) 
A ( 3610 ) 
A ( 361 1 ) 
A ( 3622 ) 
A ( 3736 ) 
A ( 3759 ) 
A ( 3898 ) 
A ( 3909 ) 
A ( 3910 ) 
A ( 391 1 ) 
A ( 3922 ) 
A ( 4036 ) 
A (4059) 
A ( 4198 ) 
A ( 4209 ) 
A ( 42 10 ) 
A (4211) 
A( 4222) 
A ( 4336 ) 
A ( 4359 ) 
A (4498) 
A ( 4509 ) 
A ( 4510 ) 
A (451 1 ) 
A ( 4522 ) 
A (4636) 
A ( 465 9 ) 
A (4777) 
A ( 4788 ) 
A { 4789 ) 
A ( 4790 ) 
A (4857) 
A <486 8 ) 
A (4869) 



= 4. 

= 2 . 

= 0 . 

^ 0 . 

= 3. 

= 4. 

= - 18 . 

= 8 . 

= 3, 

= 0 . 

= 0 . 

= 3. 

= 4. 

= - 18 . 

= 8 . 

= 3. 

= 0 . 

= 0 . 

= 3. 

= 4. 

= - 18 . 

= 8 . 

= 3. 

= 0 . 

= 0 . 

= 3. 

= 4. 

= - 18 . 

= 8 . 

= 3. 

= 0 . 

= 0 . 

= 3. 

= 4. 

= - 18 . 

= 8 . 

= 3. 

= 0 . 

= 0 . 

= 3. 

= 4. 

= - 18 , 

= 8 . 

= 3. 

= 0 . 

= 0 . 

= 3. 

= 4. 

= - 18 . 

= 8 . 

= 1 . 

= 4. 

= - 10 . 



INPUT THE BOUNDARY CONDITIONS 



C ( 1 ) = 
C ( 2 ) 1 
C ( 3 ) = 
C ( 4 ) = 
C ( 5 ) : 
C C 6 ) : 

cm = 

C ( 8 ) : 
C ( 9 ) = 
C ( 10) 
C { 1 1 ) 
C ( 12 ) 
C( 13) 
C ( 25 ) 
C ( 37 ) 
C ( 45 ) 



—46 * 
-32. 
-28. 
-24. 
- 20 . 
-16. 
- 12 . 
-2.5 
- 2 . 

* -1.5 
= - 1 . 

= -.5 
= - 10 . 
= - 10 . 

• -io. 
= -16. 
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ooooo ooooo 



180 



210 



220 



230 



240 



C ( 49 ) = -50. 

C ( 56) = -6. 

C ( 57) = -8. 

C (60) = -10. 

C ( 67 ) » -6.25 
C ( 68) = -3.125 
C(7S) = -10.5 
C ( 79) = -14. 

C ( 90) = -3.5 
C ( 71 ) = -50. 

C ( 82 ) = -10. 

10 . 

-30. 
- 10 . 
- 10 . 

10 . 
- 10 . 
- 10 . 
- 10 . 
-19. 

- 3 . 

-7. 

- 6 . 

-5. 

-9. 

-2.5 
-2. 

-1.5 
-1. 

-0.5 
= 1,201 
C ( I ) 



SOLVE THE TWO HUNDRED AND GNE EQUATIONS 
SIMULTANEOUSLY TO OBTAIN A SOLUTION FOR THE 
STREAM FUNCTION AT DISTINCT POINTS IN THE FLOW. 

CALL GELB(B, A, 201, 1,12,12, . 001 , I E R ) 

WRITE (6,900) 

DO 210 I = 1,201 
D ( 1 ) = B ( I ,1) 

CONTINUE 

WRITE THE SOLUTION TO THE STREAM FUNCTION AT 
EACH DISTINCT POINT ALLOWING FOR VARIABLE MESH 
CONSIDERATIONS 

DO 220 I = 1,37,12 
J = I 

K = J + 11 

WRITE ( 6 ,400 ) (D(L),L=J,K) 

WRITE (6,300) 

CONTINUE 

DO 230 I = 49,71,11 
J = I 

K = J + 10 

WRITE(6» 500 5 (D(L),L = J,K) 

CONTINUE 
WRITE (6,300) 

DO 240 I = 82,94,12 
J = I 

K = J + 11 

WRITE (6 ,400 ) (D(L) ,L=J»K) 

WRITE (6,300) 

CONTINUE 

DO 250 I = 106,190,12 
J = I 

K = J + 11 

WRITE (6,400) (D(L),L=J,K> 

WRITE (6,900) 



C ( 106 ) 

C ( 1 18 ) 

C ( 130 ) 

C ( 142 ) 

C( 154) 

C( 166) 

C ( 178 ) 

C ( 190 ) 

C( 191 ) 

C ( 192 ) 

C ( 193 ) 

C( 194) 

C( 195) 

C ( 196 ) 

C ( 197 ) 

C( 198) 

C( 199) 

C ( 200 ) 

C ( 201 ) 

DO 180 
B (I ,1 ) = 
CONTINUE 



I 
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250 CONTINUE 
STOP 
END 



A sample output follows on the subsequent page. For 
this case, the flow velocity was evenly distributed in the 
canard augmentor and was constrained to be twice the 
velocity of freestream conditions. 

Boundary conditions are not included in the output. 

The decimal point of each stream function value represents 
the approximate relative position of each solution point in 
the domain depicted in Figure 2. 
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